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Abstract. We report on recent results about entanglement in the context of particle mixing 
and oscillations. We study in detail single-particle entanglement arising in two-flavor neutrino 
mixing. The analysis is performed first in the context of Quantum Mechanics, and then for the 
case of Quantum Field Theory. 



1. Introduction 

Entanglement has been widely investigated in a number of physical systems, ranging from 
condensed matter to atomic physics, and quantum optics [Ij. Also in the context of particle 
physics, the role of entanglement has been considered, see for instance Refs. [2]. 

In this paper, we consider the entanglement associated to neutrino mixing and oscillations. A 
detailed study of such a topic has been performed recently in the context of quantum mechanics 
O Uj. Here we review the main results of these studies in the simplest case of two flavors. We 
show that these results suggest a simple extension of the analysis to the relativistic domain, thus 
providing the physical basis and the mathematical tools for the quantification of entanglement 
in quantum field theory. 

The phenomenon of particle mixing, associated with a mismatch between flavor and mass of 
the particle, appears in several instances: quarks, neutrinos, and the neutral ii'-meson system 
[SI E]. Particle mixing is at the basis of important effects as neutrino oscillations and CP 
violation [7]. Flavor mixing for the case of three generations is described by the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) in the lepton instance [8.) 9]. 

In the following we consider only the simplest case of two flavors. In such a case, the PMNS 
matrix reduces to the 2x2 rotation Pontecorvo matrix U(0), 

/ cose sine\ 
^ ^ \ — sm cos (y J 

which connects the neutrino states with definite flavor with those with definite masses: 

\iAf)) = u{e) |z£(™)) (2) 

where \iAf^) = {\u,),\iy^)f and \iy^"'^) = {\iyi),\u2)f. 



From Eq. ([2]), we see that each flavor state is given by a superposition of mass eigenstates, 
i.e. I fa) = Uai\i'i) + Ua2\v2)- Let us recall that both and are orthonormal, i.e. 

{l^a\l'l3) = 6a,l3 and {l^ill^j) = 5ij. 

We now establish the following correspondence with two-qubit states: 

= |l)l|0)2 = |10), |Z.2) = |0)l|l)2 = |01), (3) 

where \)i denotes states in the Hilbert space for neutrinos with mass m^. Thus, the occupation 
number allows to interpret the flavor states as constituted by entangled superpositions of the 
mass eigenstates. Quantum entanglement emerges as a direct consequence of the superposition 
principle. It is important to remark that the Fock space associated with the neutrino mass 
eigenstates is physically well defined. In fact, at least in principle, the mass eigenstates can be 
produced or detected in experiments performing extremely precise kinematical measurements 
|10] . In this framework, as discussed in Ref. [3j, the quantum mechanical state ([2]) is entangled 
in the field modes, although being a single-particle state. 

Mode entanglement defined for single-photon states of the radiation field or associated 
with systems of identical particles has been discussed in Ref. [11]. The concept of mode 
entanglement in single-particle states has been widely discussed and is by now well established 
[T2] . Successful experimental realizations using single-photon states have been reported as 
well |13) . Moreover, remarkably, the nonlocality of single-photon states has been experimentally 
demonstrated verifying a theoretical prediction [15]. Furthermore, the existing schemes to 
probe nonlocality in single-particle states have been generalized to include massive particles of 
arbitrary type [16j. 

In the dynamical regime, flavor mixing (and neutrino mass differences) generates the 
phenomenon of neutrino oscillations. The mass eigenstates {vj) have definite masses rrij and 
definite energies LVj. Their propagation can be described by plane wave solutions of the form 
= e~^^^*\uj). The time evolution of the flavor neutrino states Eq.([2]) is given by: 



= U(t)|z.(^)) , U{t) = V{e) Uo(t) IJ{0)-' , (4) 

where \iA^^) are the flavor states at t = 0, Uo(i) = diagie'^'^^^ e'^"^^*), and U(t = 0) = I. 
At time t, the probability associated with the transition i/q, — t- i/^ is 

Pu^-^u,{t) = K^/3|^«(t))P = |Ua/3(t)|', (5) 

where a, /3 = e, /j, . The explicit form for the transition probabilities in the two flavor case is: 

= l-sin2(2e)sin2(^^^^t), (6) 
P.,^.,(t) = sin2(20)sin2(^^^t). (7) 



Flavor neutrino states are well defined in the context of Quantum Field Theory (QFT), where 
they are obtained as eigenstates of the flavor neutrino charges |17|, I18j. In the relativistic limit, 
the exact QFT flavor states reduce to the usual Pontecorvo flavor states Eq.Q: flavor modes 
are thus legitimate and physically well-defined individual entities and mode entanglement can 
be defined and studied in analogy with the static case of Ref.[3j. We can thus establish the 
following correspondence with two-qubit states: 

= |l)e|0)^, = |0)e|l>/.. (8) 



States \0)a and |1)q, correspond, respectively, to the absence and the presence of a neutrino 
in mode a. Entanglement is thus established among flavor modes, in a single-particle setting. 
Eq. ^ can then be recast as 

Kit)) = U„e(t)|l)e|0);, + Ua;.(t)|0)e|l)^ , (9) 

with the normalization condition X)/? = 1 (a, /3 = e, /x). The time-evolved states 

\ii^-^\t)) are entangled superpositions of the two flavor eigenstates with time-dependent 
coefficients. Thus, flavor oscillations can be related to bipartite (flavor) entanglement of single- 
particle states [4j. 

2. Entanglement in neutrino mixing — Quantum Mechanics 

As discussed in the Introduction, the flavor neutrino state at a given time, say \i^e(t)) for 
definiteness, can be regarded as an entangled state either in terms of the mass eigenstates 
or in terms of the flavor eigenstates (at a fixed time). In the first instance, which was studied in 
detail for the multipartite case in Ref.[3], we have a static entanglement, in the sense that the 
result of the entanglement measures on the state |fe(i)) do not depend on time. In the second 
case, considered for the general three flavor case in Ref.^, the entanglement varies with time 
as it is related to the oscillations of fiavor(s). 

In this Section, we discuss these two forms of entanglement in the simple case of two fiavors, 
by means of entropy measures first, and then using the characterization of entanglement in terms 
of quantum uncertainties. The second approach turns out to be suitable for the generalization 
of our discussion to the case of mixing among relativistic quantum fields - see next Section. 

2.1. Neutrino entanglement via linear entropy 

In terms of the mass eigenstates, the electron neutrino state at time t reads: 

\u^{t)) = e-''^'^ cos9\iJi) + e-'"^^* sm9\u2) , (10) 

where {ui) are interpreted as the qubits, see Eq.Q. 

Following the usual procedure, we construct the density operator = \iya{t)){va{t)\ 

corresponding to the pure state |fQ,(i)). Then we consider the density matrix p-"^ = Trjlp^"^] 
reduced with respect to the index j. For the specific case of the state Eq. pU]) . we have 

p(^) = \ue{t))K{t)\ and 

= rr2[p(")] =cos2^|l)ii(l| + sin2 0|O)ii(O| (11) 
= Tn =cos2e 10)22(01 + sin2^ 11)2 2(11 (12) 
It is then easy to calculate the corresponding linear entropies, which turn out to be equal: 

Si-''\pe) = 2 (1 - Tr,[{p['^f]) = sm\2e), (13) 
5f ^)(pe) = 2 (1 - rr2[(pf ^)2]) = sm\29), (14) 

Similar results are found for the muon neutrino state. Note that the above results are 
particular cases of the more general ones obtained for the three fiavor neutrino states in Ref . [3j , 
where it was found that such states can be classified as generalized W states. In the present 
(two-flavor) case, the form of the entangled state is simply that of a Bell state. 

Eqs. ()13p - ()14p express the fact that flavor neutrino states at any time can be regarded as 
entangled superpositions of the mass qubits where the entanglement is a function of the 
mixing angle only. 



Figure 1. (Color online) Linear entropy S^' (pe) (full) as a function of the scaled time 
T = ■ The mixing angle 9 is fixed at the experimental value sin^ 9 = 0.314. The transition 

probabilities Py^^y^ (dashed) and P^^^^^ (dot-dashed) are reported as well for comparison. 



Let us now turn to the dynamic entanglement arising in connection with flavor oscillations. 
To this aim, we rewrite the electron neutrino state \U(.{t)) as 

\Ve{t)) = \Jee{t)\Ue) +\^et,{t)\v^) , (15) 

where |i^e)j are the flavor neutrino states at time i = and are now taken as the relevant 
qubits (see Eq.([8])). By proceeding in a similar way as done for the static case, we arrive at the 
following expression for the linear entropies associated to the above state: 

St'\p,) = S^l-'^\p,) = 4|Uee(t)nUe,(t)P 

= 4|Uee(t)|'(l-|Uee(t)|') (16) 

Eq. (|16p establishes that the linear entropy of the reduced state is equal to the product of the 
two- flavor transition probabilities given in Eqs.([5])-(I7]). It is remarkable that simple expressions 
similar to those of Eq. ()16p hold also for the three flavor case |4j . 

Note also that, for any reduced state p of a two-level system one has that 5^ = 2[1 — Tr(/j^)] = 
ADetp = 4Ai(l — Ai), where Ai is one of the two non-negative eigenvalues of p, and the 
relation Ai + A2 = 1 has been exploited. Comparing with Eq. ()16p . one sees that the transition 
probabilities coincide with the eigenvalues of the reduced state density matrix. 

In Fig. [T]we show the behavior oiS^^'^\pe) as a function of the scaled, dimensionless time 
T = ?^^2 ■ In the same figure, we also report the behavior of the transition probabilities 
Pi/e^ue and Py^^y^. The plots have a clear physical interpretation. At time T = 0, the 
entanglement is zero, the global state of the system is factorized, and the two flavors are not 
mixed. For T > 0, flavors start to oscillate and the entanglement is maximal at largest mixing: 
Pu^^ve = -fVe^i^M = 0-^' ^'^'^ minimum at T = vr. 

2.2. Neutrino entanglement via uncertainties 

An alternative characterization of entanglement can be given in terms of quantum uncertainties 
I22j . In order to apply this formalism to the case of neutrino mixing and oscillations we 
introduce the (fermionic) annihilation operator Ui for a neutrino with mass m^, with anti- 
commutator {aj,aj} = 5ij. We then define neutrino states with definite masses as: 

= aj|0)^, i = l,2 (17) 



where |0)m = |0)i (E> |0)2 is the vacuum for the mass eigenstates. 

Next we define the flavor annihilation operators by means of the following relations: 

ae{t) = cos6ai{t) +sm6a2{t) (18) 
Oi^{t) = — sinO a2{t) + cos9ai{t) (19) 

where ai{t) = e^^^^ai, with i = 1,2. 
The flavor states are given by: 

\u,{t)) = ai{t)\Q)m, a = e,iJ. (20) 

We use in the following the notation ji/g-) = \t^a{t = 0)). 
The hamiltonian for the system is 



H 



wiajai + 002^2^2 (21) 



where we used the relations coee = cos^ 9 + U02 sin^ 6 , lo^^ = loi sin^ 9 + L02 cos^ 6 , iOei_i = 
{u}2 — ^1) sin COS 9 . 

In the mass basis, we can introduce the SU{2) operators (the superscript m stands for mass): 

= ajaa , J" = a^ai, (22) 
= \ («i«i - «2«2) (23) 



C 



^ {a\ai + a\a2^ (24) 



We also have 



Ni = C + (25) 
N2 = C - (26) 
H = u)iNi + U2N2 (27) 

The static entanglement of the electron neutrino state \ve{t)) defined in Eq. ()20p . is 
characterized, in the present formalism, by the variances associated with the numbers A'j, relative 
to the mass qubits. Thus we have: 

^Ni{Ve) = {Ue{t)\Nf\Ve{t)) " (l^e (t) | iVi | ^^e (*))' 

= ^sin2(2e), i = l,2. (28) 

This result differs by a factor 4 from that obtained by means of the linear entropy, Eqs. (|13|) - (|14p . 

In order to discuss the dynamical entanglement of the state \V(,{t))^ we need to introduce 
flavor oscillations, which can be seen both in terms of overlaps of states at different times: 

Pu^^uAt) = \{^e\yem'' (29) 

P..^u,{t) = \{y^v,{t))\'' (30) 



with P^^^^^{t) + Py^-^u^{t) = 1) or equivalently in terms of expectation values of number 
operators at time t: 



Pu,-,uM = {l^e\Neit)\l^e) 



(31) 
(32) 
(33) 



The explicit expressions of the transition probabilities are given in Eqs.([6]),(I7]). 
In the flavor basis, we can again introduce SU{2) operators (at time t): 

Jlit) = al{t)a^{t) , j[{t) = al(t)ae{t) 

4{i) = \ (al(t)ae(t) - «|.(t)«M(*)) 



where the superscript / stands for flavor. The flavor number operators are: 

NS) = C + 4{t) 
N,{t) = C- 4{t) 

Note that the Hamiltonian can be written as: 

H = LOeeNeit) + iOf,^N^{t) + iOe^,{ji{t) + jZ(t)) 



(34) 
(35) 

(36) 



(37) 
(38) 



(39) 



Flavor entanglement is given by the variances of the above flavor numbers. Consider first the 
the variances of the SU (2) operators in the flavor basis, which give: 



AJl(z.e)(t) 

AJ2(j^e)(i) 
AJ-i{Ue){t) 



- 1 — sm (46') sm I 

^ - sin2(26l) sin2[(w2 - wi)i] 



sin2(26')sin' 




2 f UJ2 — iOl 



1 — sm [10) sm ] 



(40) 

(41) 

(42) 
(43) 



When we restrict to a given flavor (e.g. electron neutrinos), we find that the entanglement is 
given by 



ANe{Vem = {u,{t)\Nl{t)\Ve{t)) - {Ve{t)\Ne{t)\Ve{t)f 



P. 



MPu^^uAt) 



(44) 
(45) 



with the same result for A.N^{ve){t). The above result coincides (again up to a factor 4) with 
the one obtained in Eq. ()16p by means of the linear entropy. 
Analogous results are easily obtained for the state |f^(t)). 



3. Neutrino mixing in Quantum Field Theory 

We now look for an extension of the above discussion to a relativistic context. To this 
aim, we need to consider neutrino mixing at level of fields pT]. For two flavors, the mixing 
transformations are 

Vei^) = cos^z/i(x) +sm6i>2{x) (46) 
v^ix) = — sin6'z^i(x) + cos 6* 2^2(2;) , (47) 

where Veix) and f/i(x) are the Dirac neutrino fields with definite flavors. Here, z^i(x) and 2^2(2^) 
are the free neutrino fields with definite masses mi and m2, respectively. For the purpose of 
discussing the quantization of fiavor fields and the phenomenon of neutrino oscillations, it is 
sufficient to consider only the free Lagrangian term for these two free fields: 

Cu{x) = vn,{x) (i^- Mi) v^{x) , (48) 

where v'^ = (1^1,2^2) and = diag{mi,m2). Ci,{x) is invariant under global U{1) phase 
transformations of the type i^mi^) — ^^°^^m{x)- This implies the conservation of the Noether 
charge Q = J I^{x)d'^x (with I^ix) = h'm{x)^^Vrn{x)) which is indeed the total charge of the 
system, i.e. the total lepton number of neutrinos. 

Consider now the global SU{2) transformation |19) : 

u'^{x) = e'''r^^iym{x) j = 1,2,3. (49) 

with aj real constants, tj = crj/2 with aj being the Pauli matrices. 

Cu is not invariant under the transformations ()49p since mi 7^ m2. By use of the equations 
of motion, we obtain 

dCu = iajUm{x) Tj, i^mix) = -ajd^J'^ ,j{x) , (50) 
where the currents are: 

= ^m{x) r Tj Vm{x) , j = 1, 2, 3. (51) 

The related charges 

Qmj(0 = 1 d3xj0_^.(x), (52) 

satisfy the sii(2) algebra: 

[Qm,i(i);Qm,i(*)] = ■i£ijkQm,k{t)- (53) 

The Casimir operator is proportional to the total (conserved) charge: Qm,o = \Qu and also 
Qm.3 is conserved, due to the fact that is diagonal. This implies that the charges for ui and 
V2 are separately conserved. The U{\) Noether charges associated with vi and V2 can be then 
expressed as 

Ql = \q + Qm,3 ; Q2 = ^Qu - Qm,3 ■ (54) 

Qi = (f^vl{x) Vi{x) , (55) 



with Q total (conserved) charge and i = 1, 2. 

Let us now consider the Lagrangian Cy{x) written in the flavor basis 

= Uf{x) {i^- M^) Uf{x) , (56) 

where = (ue, v.) and = ( "^"^ "^^^^^ V 

The variation of the Lagrangian (j56p under the SU{2) transformation: 



is given by 



where 



Again, the charges 



u'j{x) = e'''^-^^uf{x) j = 1,2,3, (57) 



6Cuix) = iajVf{x) [Tj, Ml,] Vf{x) = -ajd^J'^-{x) , (58) 



J^^ix) = Uf{x) 7" Tj Uf{x) , i = 1, 2, 3 . (59) 



QfAt) = J d'^Jh(^) (60) 



close the su{2) algebra, however, because of the off-diagonal (mixing) terms in M^, Qf^t) 
is time dependent. This implies an exchange of charge between and v^, resulting in the 
phenomenon of neutrino oscillations. The (time dependent) flavor charges for mixed fields are 
then defined as [19]: 

Qe{t) = \Q + Qf,z{t) , Q^.{t) = \Q-Qf,^{t) , (61) 
Q„{t) = J d^xi^lix) v^x) , (62) 

where a = e, fj, and Qe{t) + Q^{t) = Q . 

3.1. Flavor states for mixed neutrinos 

Till now our considerations have been essentially classical. We now quantize the fields with 
definite masses as usual (see Appendix), and consider the eigenstates of the above defined 
charges. 

The normal ordered charge operators for free neutrinos vi, V2 are: 

■.Q,:^ /rf'x : v\{x) v,{x) := J d^k [a'^^^al^, - /Sl/^k,) , (63) 

where i = 1,2 and : .. : denotes normal ordering with respect to the vacuum |0)i^2- The neutrino 
states with definite masses defined as 



k,i/ 



a;;^JO)i,2, i = l,2, (64) 



are clearly eigenstates of Qi and Q2, which can be identified with the lepton charges of neutrinos 
in the absence of mixing. 



The situation is more delicate when mixing is present. In such a case, the flavor neutrino 
states have to be defined as the eigenstates of the flavor charges Qo-(t) (at a given time). The 
relation between the flavor charges in the presence of mixing and those in the absence of mixing 
is: 

Qe{t) = cos'^d Qi+sin'^0 Q2 + sin^cosS / d^x\ul{x)u2{x) + i'2{x)i'i{x)\ , (65) 
Qf,{t) = sm"^ 9 Qi+cos^e Q2- smOcose J d^x i^l{x)i^2{x) + i'l{x)ui{x) . (66) 

Notice that the last term in these expressions is proportional to the charge Qm,i defined above 
(cf. Eq. (j52p l. The presence of such a term forbids the construction of eigenstates of the Qo-(t) in 
the Hilbert space ?^i.2- One then is led [17J to define another Hilbert space, Tie.fj, for the flavor 
neutrino fields. The flavor vacuum state |0)e,/x and the mass vacuum state |0)i^2 are orthogonal 
to each other (see Appendix). 

The normal ordered flavor charge operators for mixed neutrinos are then written as 



:: Q„{t) :: = j d^x :: i/t(x) u„{x) :: 

r 

where o" = e,/i, and :: ... :: denotes normal ordering with respect to |0)e,^. Thus, the flavor 
charges are diagonal in the flavor annihilation/creation operators constructed by means of the 
mixing generator Gg defined in the Appendix. The definition of the normal ordering :: ... :: for 
any operator A, is the usual one: 

■.■.A:: = A - ,^f,{0\A\0)e,^,. (68) 

Note that :: Q„{t) :: = \t) : Qj : Gg{t), with (aj) = (e, 1), (^, 2), and 

:: Q :: = :: Qe{t) :: + :: g^(t) :: = : : + : Q2 : = : Q : . (69) 

We define the flavor states as eigenstates of the flavor charges Qa at a reference time t = 0: 

l'^k,-)=«i(0)|0(0))e,^, a = e,M (70) 
and similar ones for antiparticles. We have 

QeiO) :: Wl,,) = Wl^e) , Q,iO) Wl,,) = Wl,,) (71) 



" Qe{0) :: l^^'k.^) = = " Qm(0) " Wl,e) , ■■ QAO) :: |0)e,;. = 0. (72) 
The explicit form of the flavor states If^e) and jf^^) at time t = is given in the Appendix. 

3.2. Oscillation formulas in Q FT 

Flavor oscillation formulas are derived by computing, in the Heisenberg representation, the 
expectation value of the flavor charge operators on the flavor state. We have 



{<.(i),<e(0)}f + |{/3:UW'«k!e(o)}f (73) 



and 



e,;.(0| :: Qe{t) :: \0)e,^ = e,;.(0|Q^(t)|0) 

The oscillation formulas are |18j : 

2 / ^k,2 — ^k, 



Qu^^uAt) = l-sin^(2e) iC^klsin 



t +|yk|^sin^ 



2 ■ 2 f ^k,2 + Wfc^i 



/-|k /.\ • 2/r,/)N r irr |2 ■ 2 f ^k,2 — ^k,l 

Qu,-^u^[t) = sm (20)^1 C/kl sm I ' ^ ' ' 

The charge conservation is ensured at any time: 



t] +|Fkl sin^ 



2 ■ 2 ( ^k,2 + l^fc,l 



(74) 

(75) 
(76) 

(77) 



The differences with respect to the Pontecorvo formulas are: the energy dependence of the 
amplitudes, and the additional oscillating term. In the relativistic limit: |k| ^> y/mim2, we have 



1 and iVkl 



and the traditional formulas are recovered. 



4. Entanglement in neutrino mixing — Quantum Field Theory 

Following what done in the QM case, we now calculate the entanglement associated to an electron 
neutrino state at time t, by means of the variances of the above discussed charge operators. 

Let us start with the Noether charges Q^- , which are expected to characterize the amount of 
static entanglement present in the states Eq. ()70p . We obtain: 



AQi{Ue)it) 



]sm\2e), i = l,2. 



r \2 

el 



(78) 



in perfect agreement with the quantum mechanical result Eq. ()28p . 

Next we consider dynamic entanglement, which is described by the variances of the flavor 
charges. We have: 



r \2 
el 



k,el 



H 2 



k,e/ 



^ / d^p (a;yi)a^,,(t) - /3^;,e(t)/3ip,e(i) 

= (^'k,el<eW«k,e(i)l^'k,e) + k, J/^^e (0/3^k,e (*) k'k.e) 

We now consider the third term in Eq. (|79p . We have: 

{^\Mie^)0^\.At)f^Xe{tW-Ke^^^^^^ 

{«lc,eW,<e(0)}|'e,M(0|/3:U(*)/3-k,eW 

- {«L,e(0),/3:k,eW}f eA^Wteii)'AAt) 



(79) 



+ 



- {<(0),/3:l,.(t)} {«k,e(0),aL!eW} eA^Wl^e^Wit) 



(80) 



Figure 2. QM vs. QFT flavor entanglement for |fe(i)) as a function of the scaled time T 
with 6 fixed at the value sin^ 9 = 0.314. 
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Explicit calculation of the above quantity shows that: 



(^'k,eKe(i)«k,e(t)/3\e(0/3:k,e(i)k'k,e) 



(^'k,el/3:l,e(0/3:k,e(i)k'k,e) 



(81) 



SO that we have 



(82) 



which formally resembles the quantum mechanical result Eq. ()44p . The differences are now due 
to the presence of the fiavor condensate, which affect the oscillation formulas (see Eqs. ()75p . ()76p ). 
In Figl^l fiavor entanglement formula is plotted in the QFT case against the corresponding QM 
case. 

5. Conclusions 

On the basis of recent results, we have discussed some aspect of entanglement in the phenomenon 
of neutrino mixing and oscillations. In the simple case of two fiavor mixing, we have shown how 
to generalize previous results obtained in the context of quantum mechanics, to the case of 
quantum field theory. The difference between the QFT and the QM cases are related to the 
condensate vacuum structure associated to neutrino mixing in QFT. 

Our study provides a simple, exactly solvable example for a possible extension of entanglement 
to quantum field theory. Apart from the differences with QM due to the fiavor vacuum 
contributions, the QFT result is interesting from a more conceptual point of view. Indeed, 
both the static and the dynamical entanglement arise in connection with unitarily inequivalent 
representations: in the case of the static entanglement, the fiavor Hilbert space at time t to 
which the entangled state {I'fjlt)) belongs, is unitarily inequivalent to the Hilbert space for the 
qubit states [17]; on the other hand, in the case of dynamical entanglement, where the 
qubits are taken to be the fiavor states at time t = 0, the inequivalence is among the flavor 
Hilbert space at different times [23]. In the first case, the relevant orthogonality relation is 
limv^oom(0|0(t))/ = 0, in the second limy^oo /(0(t')|0(t))/ = 0, with t ^ t' . 

Since the inequivalent representations are associated with a non-trivial condensate vacuum 
structure, the above conjecture suggests that, in the context of QFT, many interpretational 
issues connected with entanglement could be revisited in this new light. 



Appendix A. QFT formalism for mixed fields 

The fields I'dx) and t'^(x) are defined tlirougli the mixing relations ()47p . in terms of the free 
fields 1^1 (x) and 1^2(2;) which are expanded as 



1,2 



(A.l) 



with aj;^,(t) = aj;^i e-*'^".'*, /3;;|.(t) = /3;;|. e*"^".'*, and ioj,^, = ^k^ + ml The operator ^ and /3|; ., 
i = l,2, r = l,2 are the annihilator operators for the vacuum state |0)m = 10)1(8)10)2: «k j|0)m = 
/^kJO)m = 0. The anticommutation relations are: jt'f (x), = 6^{'x. — y)dai36ij, 

with a,f3 = I, ..A, and {ak^j,aq|j| = ^kq'^rs'Jjj; {/^k,,, /?qj} = '^kq'Jrs'^ij, with ij = 1,2. All 
other anticommutators vanish. The orthonormality and completeness relations are given by 



ul^^ut , = vl^vt , = 6rs, ulW-k.i = v^juIj^ = 0, and Er(«k,j<!i + vLuAd = I. 



The generator of the mixing transformations is given by [17j: 



Gg{t) = exp 



d x I i/l{x)i>2{x) — z^2(a^)'^i(a^) 



(A.2) 



so that 



i^^ix) = Gg'it) ut{x) Ge{t) ; {a,i) = (e, 1), (^,2) 



(A.3) 

At finite volume, this is a unitary operator, Gg^{t) = G-0{t) = G\{t), preserving the 
canonical anticommutation relations. The generator Gg^{t) maps the Hilbert space for free 
fields Hi^2 to the Hilbert space for mixed fields 7ie,fj.- Gg^{t) : Hi^2 ^ T^e,fj.- In particular, the 
flavor vacuum is given by |0(t))e,/x = Gg^{t) |0)i,2 at finite volume V. We denote by |0)e,^ the 
fiavor vacuum at t = 0. In the infinite volume limit, the flavor and the mass vacua are unitarily 
inequivalent [T7]. The explicit expression for |0)e,;i at time t = in the reference frame for which 
k=(0,0, |k|)is 



|0) 



n[(l-sin2e|Fk|')-e"sine cos 9 iV^lia^^^^l^^ + (^K2^Xi) + (^-4) 



+ sin^ 9 mU^\ia'^^,^Xi " <2/3l,2) + \V^\' at/^X^^^I^Xi] |0)i,2 
The condensation density is given by 

e,M(0|<.«k,JO)e,;. = eA0\^i^^(3L,^\0)e,^. = sin^^ |Fk|' , ^ = 1, 2 . 

The flavor flelds are written as: 
1 



(A.5) 



^.(x,t) = <,^<^it) + ^-K^f^-iJt) ^ (cT,i) = (e,l),(^,2). (A.6) 

V ^ k,r ^ ^ 

The flavor annihilation operators are jl7j : 

"k,e(i) = cos9al^{t) + sm9 XI ['"kli'"k,2 "k,2(*) + «M'^-k,2 /3!,k,2(*) 

s 

"k,/.(i) = cos9al2{t) - sm9 XI ['"k!2^k,i "k.iW + u^,2^\i /3%^^{t) 

s 

/3L^^,{t) = cos9 (3L^^,{t) + sine XKk,2^-k,i/3ik,2W + O-k,! "kV*)" 

s 

f3\,{t) = cose /3:k,2(i) - Sine X[A,i^-k,2/3ikAW + <i^-k,2<iW]- (A.7) 



In the reference frame where k = (0, 0, |k|), we have 

ai^,{t) = cose ai^,{t) + sine {\U^\ a^^it) + \V^\ ^Xd*)) ^ 

aljt) = coseai^^it) - sine {m a^^t) - \V^\ (3X,i{t)) , 

/3:k,e(i) = cos0/3:k,i(t) + sine {\U^\ f3L^^^{t) - \V^\ a'^^^it)) , 

= cos0/3:k,2(t) - sine (|t/k| PL^^^(t) + \V^\ a'^.H)) . (A.8) 

Here = (—1)^ and 



2./a;fc,ia;fc,2(wfc,i + mi){ujk,2 + ?7i2) 



with iC/kP + l^kP = 1. 

The expUcit expressions for the flavor states If^g) and jf^^) at time t = 0, in the reference 
frame for which k = (0, 0, |k|) are 

K,e) = <e(0)|0)e,;. (A.IO) 

= cos^agi + If^kl sine - \V^\ sine ak!i<2/3-k,i] G'ki^r(^) 11 G^\e)\0)i,2, 

WLJ = «i(o)|o)e,^ (A.ll) 

= cosea22 - l^kl sine a\^^, + \V^\ sine <i«L!2/31,2] G^^M^) U Gp'(^)|0)i,2 , 

P^k 

where G{e,t) = Hp ns=i G'p,s(e, t). In these states a multiparticle component is present, 
disappearing in the relativistic Umit |k| ^ ^•m\m2 : in this hmit, since |C/k| — ^ 1 and 
|Vkp — 7- 0, the (quantum-mechanical) Pontecorvo states are recovered. 
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